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ABSTRACT 

The baryon acoustic oscillation (BAO) feature in the clustering of matter in the universe serves as 
a robust standard ruler and hence can be used to map the expansion history of the universe. We 
use high force resolution simulations to analyze the effects of galaxy bias on the measurements of 
the BAO signal. We apply a variety of Halo Occupation Distributions (HODs) and produce biased 
mass tracers to mimic different galaxy populations. We investigate whether galaxy bias changes the 
no n-linear sh i fts on the acoustic scale relative to the underlying dark matter distribution presented 
bv iSeo et al.l ()2009h . For the less biased HOD models (b < 3), we do not detect any shift in the 
acoustic scale relative to the no-bias case, typically 0.10% ± 0.10%. However, the most biased HOD 
models (6 > 3) show a shift at moderate significance (0.79% ± 0.31% for the most extreme case). 
We test the one-step reconstruction technique introduced by lEisenstein et al.1 ((2007) in the case of 
realistic galaxy bias and shot noise. The reconstruction scheme increases the correlation between the 
initial and final (z = 1) density fields achieving an equivalent level of correlation at nearly twice the 
wavenumber after reconstruction. Reconstruction reduces the shifts and errors on the shifts. Wc find 
that after reconstruction the shifts from the galaxy cases and the dark matter case are consistent with 
each other and with no shift. The la systematic errors on the distance measurements inferred from 
our BAO measurements with various HODs after reconstruction are about 0.07% — 0.15%. 

Subject headings: distance scale — cosmological parameters — large-scale structure of universe — 
baryon acoustic oscillations — methods: N-body simulations — cosmology: theory 



1. INTRODUCTION 

Baryon acoustic oscillations (BAOs) are the im- 
pressions of acoustic waves in the hot plasma of the 
early universe onto the matter distribution. The 
acoustic waves were imprinted onto the large-scale 
baryon distribution once the photons and baryons 
decoupled after recombination with a characteristic 
wavelength of about 150 Mpc. This characteristic or 
acoustic scale is the distance traveled by the sound 
waves before recombination and is known as the "sound 
horizon scale" , which has also been measured in the 
Cosmic Microwave Background (CMB) to be about 
one degree. Since galaxies form in overdense regions 
of the universe, they are preferentially separated by 
this acoustic scale. The acoustic scale can therefore 
be used as an excellent standard ruler. More recently, 
the BAO scale has been measured in large scale sur- 
veys as a peak in the correlation function £(r) or 
as a damped harmonic series in the power spectrum 
P(k) (lEisenstein et al.ll2005HCole et a 
Tegmark et al.l l2006t IPercival et al] 
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there has been growing interest in using the BAOs as a 
powerful probe of dark energy. BAO measurements have 
become a key component of many dark energy surveys, 
both current and future. 

The systematic errors on measurements of the acous- 
tic scale are dominated by non-linear structure growth at 
small scales and low redshift and by redshift-space dis- 
tortions. Non-linear effects broaden and shift the peak 
in the correlation function. Hence, it is harder to cen- 
troid the peak and obtain precise and accurate mea- 
surements of the location of the peak. Redshift dis- 
tortions further dampen the BAO signal along the line 
of sight. The largest dark energy surveys plan to mea- 
sure the acoustic scale to the 0.1% — 0.2% level. While 
most of the growth on BAO scales (~ 150 Mpc) is still 
linear, to achieve percent and sub-percent level accu- 
racy, it becomes imperative to understand and model 
the shift and systematic errors caused by non-linear 
effects. There has been an increasing effort in un- 
derstanding the effects of non-linear growth and non- 
linear gravitational flows in order to model them for 
surveys. Most of the modeling has been done using 
dark matter simulations that do not inco rporate the ef- 
fects o f observing galaxie s as m as s tracers ( Mciksin 'et all 
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populate their dark matter ha los with galaxy populations 
(also see ISanchez et all (|2008l) ). However, their high res- 
olution simulation covered a volume of 2.41ft," 3 Gpc 3 . 
The precision of BAO measurements is proportional to 
the volume used. In order to simulate real surveys we 
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want to use several tens of /i~ 3 Gpc 3 . In our simula- 
tions, we use a variety of galaxy prescriptions but with 
a volume almost 20 times larger simulation, and explore 
the effects of reconstruction on the BAO signal. 

Redshift surveys observe galaxies and not the under- 
lying matter distribution. It is well known that galax- 
ies form in the high density regions of dark matter ha- 
los. Thereby using galaxies as mass tracers introduces 
an inherent bias into any survey by over-weighting the 
overdens e regions and under- weighting the underdense 
regions. iPadmanabhan fe White! (|2009T l use perturba- 
tion theory to show that we should expect the shift on 
the acoustic scale to be less than a percent even with 
the presence of galaxy bias. In order to measure this 
effect precisely, we need to run simulations with t he ef- 
fects of galaxy bias. While iSeo et all (|2008l 12009^ show 
the shift in the acoustic scale due to nonlinear structure 
growth and redshift distortions, here wc include the effect 
of galaxy bias in our highest resolution simulation set. In 
this work, we use a set of high resolution simulations with 
various Halo Occupation Distributions (HODs) to model 
the galaxy p opulations expected to b e measured by vari- 
ous surveys (jZheng et al.ll2007l I2009T ) . We investigate the 
effects of galaxy bias on the acoustic scale for each case. 
We compare our resu lts with galaxy bias to the results 
from lSeo et al.l (|2009D . who used all dark matter particles 
t o trace the matte r densi ty field (hence, no galaxy bias). 

lEisenstein et al.1 (j2007j ) presents a simple, one-step re- 
construction technique to account for the effects of large- 
scale gravitational bulk flows that degrade the acoustic 
signal. This method has been proved to be successful 
for dark-ma tter simulations even for non-negligible shot 
noise levels (lEisenstein et al.ll2007t ISeo et al.l l2008l [200l 
INoh et al.ll2f)09t ). iHuff et al.l 1)20071) look at the effect of 
galaxy bias on reconstruction on a smaller volume (1.1 
hr 3 Gpc 3 ). In this paper, we use a much larger volume 
of 44 h~ 3 Gpc 3 to investigate the effect of galaxy bias on 
the technique of reconstruction. 

In § [2] we describe the simulations, the Halo Occupa- 
tion Distributions, fitting methods and the reconstruc- 
tion technique. We introduce the propagator and ana- 
lyze the effects of galaxy bias and reconstruction on it in 
§[3J Next, wc compute the acoustic scale and investigate 
the effects of galaxy bias on the acoustic scale in § 2] Wc 
also scompare our results t o the Fisher Matrix estim ates 
using the code provided bv lSeo fe~ 
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to the dark matter only analysis of 
summarize our results in § [5] 



2. SIMULATIONS & METHODS 

2.1. Simulations and Code 

We generate high force res olution ,/V-body sim u lation s 
using the ABACUS code by iMetchnik fe Pintol (f2010h . 
The ABACUS code uses a new method to compute N- 
body forces for periodic boundary conditions by comput- 
ing the force on a particle from the rest of the simulation 
and an infinite sum of periodic images. The force on 
a particle is calculated by dividing the simulation into 
near field and far field. The near field force is computed 
by a direct 0(N 2 ) method with a gravitational softening 
length of 0.0488/i -1 Mpc. The far-field force is computed 
by summing the forces resulting from the multipole ex- 
pansion of the mass distribution within a grid of distant 



cells, aided by a cyclic convolution to speed the calcula- 
tion. 

We assume concordance cosm ology consistent wit h 
the WMAP5+SN+BAO results (jKomatsu et al.ll200ll : 
Q m = 0.279, n A = 0.721, h = 0.701, n b = 0.0462, 
n s = 0.96, a$ = 0.817. We used 2 nd order Lagra ngian 
pertur bation theory (2LPT) based IC code by iSirkol 
(|2005l ) to generate the initial conditions. In particular, 
we use linear theory density field at z = oo and use 2LPT 
code to generate the density field at z = 50. The initial 
(z = oo) density field is computed by back-evolving the 
z = output from CMBFAST assuming a matter domi- 
nated universe. When we compare our results to the ini- 
tial density field, we always use the initial linear theory 
density field. We evolve the density fields from z = 50 to 
z = 1 using the ABACUS code and calculate the spher- 
ically averaged power spectrum in both real space and 
redshift space. We generate 44 simulations each contain- 
ing 1024 3 particles in a l/i _1 Gpc box giving us a total 
volume of 44/i -3 Gpc 3 . The resolution of our simulations 
gives us a particle mass of 7.2 x 10 10 h~ 1 M Q . Note that 
this is the same simulation set as G1024 set presented in 
ISeo et aD (|2009l ) . Figure Q] shows the effects of structure 
evolution on the baryon acoustic oscillations in the power 
spectrum at z = 1. In this figure, the power spectrum is 
divided by the no- wiggle power spectrum P avr i(k), which 
is the matter densi ty power spectrum with out the BAO 
harmonic features (jEisenstein fc HvJll998h . We see two 
distinct harmonic acoustic peaks before the acoustic os- 
cillations get erased by non-linear structure growth on 
small scales. 

2.2. Halo Occupation Distributions (HODs) 

Perturbation theory explains how overdensities in mat- 
ter grow while dark matter simulations help us bet- 
ter understand non-linear effects. In real surveys, of 
course, we do not observe the dark matter particles 
but galaxies residing in dark matter halos. We assume 
that the galaxy population are tracers of the underly- 
ing dark matter distribution. Since galaxies form in 
the highest density regions of dark matter halos, they 
represent a biased version of the matter density. It is 
critical that we understand how biased tracers of the 
matter distribution affects the acoustic scale. To ex- 
amine the effects of galaxy populations, we generated 
various Halo Occupation Distri butions (HODs) and ap- 
plied them to our simula tions (jPeacock fc Smith] 120001 : 
iBerlind fc Weinberg! 12001 ) . The HOD models allow us 
to use various mass tracers (galaxy populations) to study 
the acoustic scale. To generate the H ODs, we use a sim- 
plc Fricnds-of-Fricnds (FoF) program (jDavis et al.lfl985h 
to identify collapsed dark matter halos with a linking 
length of bpoF = 0.16. We populate each halo with a 
central and satellite galaxies using the following formula: 

N(M) = (° HM<M cen 
gK ' \l + Poisson(M/M SQt ) if M>M cen KJ 

where M is the halo mass, M cen is the minimum mass 
needed to have a central galaxy Poisson(A//Af sat ) is 
an integer randomly chosen from a Poisson distribution 
with mean M/M sat , and M sat is the minimum mass to 
have at least one sa t ellite galaxy (iGuzik fc Seliakl 120021: 
IBerlind et all l2003t iKravtsov et al.l 120041: iZheng et al l 
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12001 . Note that P(M/M sat ) = for M < M sat . To 
test the effects of various galaxy populations, we vary 
M cen and M sat to generate 12 different HOD models 
whose properties and definitions arc detailed in Table [TJ 
By varying the mass thresholds M cen and M sat , we ob- 
tain HOD models with different number of galaxies. For 
each central galaxy threshold, M cen , we vary the satellite 
mass threshold and then slightly adjust M cen to give us 
the total number of galaxies: 2 x 10 6 , 1 x 10 6 , 3 x 10 5 and 
1 x 10 5 . The central galaxy is assigned the halo's center 
of mass position and velocity. We then randomly pick a 
corresponding number of halo particles and assign their 
positions and velocities to the satellites. 

Thus, each HOD identifies halos based on the mass 
threshold and replaces the halo by a central galaxy and 
satellite galaxies (based on the satellite threshold) and 
uses these galaxies for any analysis. We use the various 
HOD models as biased tracers of the density fields and 
explore their effects on the acoustic scale. Since we seek 
to understand the effects of various galaxy populations 
or HOD models, we use the "mass" case, where we do 
not apply any HOD, to be our b ase case. A full a nalysis 
of the mass case is presented in ISeo et ahl (|2009D . The 
different number densities corresponding to each HOD 
also allow us to explore the effect of shot noise on the 
acoustic scale. 

2.3. Power Spectra Fitting Methods 

We compute the power spectrum for e ach of the 44 sim- 
ulatio ns following the method used by ISeo et al.l ([20081 
120091) . We use \ 2 analysis to fit the spherically averaged 
power spectrum, P bs, to the template power spectrum 
P m {k): 

P ohs = B{k)P m (k/a)+A{k), (2) 

where B(k) represents scale-dependent bias and redshift 
distortions; A(k) represents anomalous power such as 
shot noise and additive terms of non-linear growth; and 
a represents the ratio of the linear acoustic scale to the 
measured acoustic scale. Therefore, a — 1 measures the 
shift in the acoustic scale i.e. the error in the distance 
inferred if one used Eq. (QJ as the template, a > 1 
implies that the measured acoustic scale is shifted to- 
wards smaller scales or larger Fourier wavenumbcrs. In 
this paper, we quote our shifts in the acoustic scale as 
a — 1 (%) in percent. Our template power spectrum 
P m {k) is constructed by modifying the linear power spec- 
trum. We use a nonlinear parameter Em that degrades 
the BAO peaks in the linear power spectrum to account 
for nonlinear effects and redshift distortions. The form 
of the template power spectrum is given by: 

Prn(k) = [Plinik) - P„wl(fc)] exp ( ~ k ^ M ^j + ^nwl(fc) 

(3) 

where Pn n (k) is the linear power spectrum and 
-Pn wi(fc) is the no- w i ggle p ower spectrum as described 
by lEisenstein k, Hul (|1998| ) . Due to the large degree of 
polynomials used to fit A(k) and B(k) (Eq. our re- 
sults are not sensitive to our value of Em- We use the 
fiducial values for Em given by the Zcl'dovich approxi- 
mation: Em 5.3/i Mpc -1 for real s pace and J -Oh M pc -1 
for redshift space. As shown by ISeo et ahl (I2008D . the 
fitting results are not sensetive over a range of Em val- 



ues (AEm = ±2/i Mpc" 1 ) due to the large degree of 
polynomials used to fit A(k) and B(k). We use a 7 th or- 
der polynomial for A(k) and a 2 nd order polynomial for 
B{k) for both real space and redshift space and fit over 
40.02/i Mpc" 1 <k< 0.35/i Mpc" 1 : 

A(k) = oq + a\k + 02k 2 + ... + a^k 7 and 

B(k) = b Q + b 1 k + b 2 k 2 . (4) 

2.4. Resampling Simulations Method 

In order to measure the scatter in a and therefore the 
scatter in the acoustic scale, we use a modified bootstrap 
method. We generate 1000 subsamples by randomly se- 
lecting M boxes without replacement out of a total of N 
simulations. We then perform a \ 2 fitting to get the best 
fit a for the individual subsamples and find the mean 
a and the scatter. The scatter in th e measured a's is 
rescaled by a factor of y/ Mj (N — M) to reflect the scat- 
ter in the mean a. For our resampling method, we have 
N = 44 and we choose M = 22 so that the rescaling factor 
is unity. To analyze the scatter in a, we fit each simula- 
tion assuming a covariance matrix of the power spectrum 
given by a Gaussian random field. This method reflects 
any non-Gaussianity in the density field in the scatter in 
the best-fit a's between different simulations. 

2.5. Reconstruction with HODs 

Large-scale gravitational forces cause large scale ve- 
locity fields between the overdensities. These velocity 
fields tend to broaden the acoustic peak and degrade the 
BAO measurement. We employ a simp l e reco nstruction 
scheme introduced by lEisenstein et al.l (|2007f) . This re- 
construction scheme is based on the Zel'dovich approx- 
imation and models the large scale velocity fields to at- 
tempt to restore the BAO information and restore infor- 
mation about the initial linear density field. 

In linear perturbation theory, we can estimate the dis- 
placements of mass particles based on the density per- 
turbations using the Zcl'dovich approximation. With re- 
construction, we estimate the bulk flows based on the 
large-scale information of the observed nonlinear density 
field and undo the large scale velocity fields to recover the 
portion of BAO that has been degraded. The Lagrangian 
particle position r can be mapped to an Eulerian particle 
position x by the displacement q: 

x = r + q. (5) 

The density in Eulerian coordinates is transformed to the 
density in Lagrangian coordinates through the Jacobian: 



po df= p dx (6) 

p_ = 1 + 8 = ^L = J -i = I (7) 

Po dx Det(cj£ + dq i /dr j ) V ' 

where po is the Lagrangian density and is the same as 
the Eulerian mean density, S is the Eulerian overdensity, 
and 5 K is the Kronccker delta function. If we expand 
the expression to linear order, it becomes: 

8 = -V • q. (8) 
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In the linear regime growing mode in CDM cosmologies, 
q is a curl-free field. Hence, we can define a scalar field 
<fi such that: 

q = -Vcf>. (9) 

In our reconstruction scheme, we are looking to undo 
large scale velocity flows and thus we smooth the small- 
scale perturbations by using a Gaussian filter. Thus, 
we enter only the large-scale information to derive the 
bulk flows. In our analysis, we use a Gaussian smoothing 
width of R = lAh" 1 Mpc to obtain the real-space galaxy 
density field <5 rcal in configuration space. Using Eq (8) 
and Eq ©, we can express <j> in terms of <5 real as follows: 



72 j,rcal 



real _ 



srrcal 



and 

d 3 k <5 rcal 1 



ik-r 



(10) 

(ii) 



where <5 rcal is the real-space smoothed density field in 
Fourier space. While we use Eq (fTT|) to solve Eq © in 
Fourier space, it can also be easily solved in configuration 
space as a central force problem. Hence, we call the 
displacement field derived from solving Eq (S) the trivial 



displacement field q 



real 
trivial 



which is given by: 



V • n real = 

y trivial 

(T roal = 
ytrivial 



real 



b 

-V0 rcal 



and therefore 



(12) 



In real space, this trivial displacement field is the desired 
displacement field. In our simulations, we displace the 
data points and the random points by the displacement 
field <?* trivial ■ By displacing both the data and random 
points, we ensure that the amplitude of the overdensities 
remains the same while removing the effects of the large 
scale velocity field. 

In galaxy redshift surveys, the density field we ob- 
serve is subject to redshift distortions. We can relate 
the redshift-spacc density field <5 rod in Fourier space to 
that in real space as follows: 



frcd 



(i + /VM 



real 



(13) 



where fi — k z /\k\ and z is the line-of-sight direction. 
We can then compute the redshift-space scalar field (j) Ted 
from Eq (II) and Eq (p): 



d 3 k S red 1 



1 



Ak-r 



(2tt) 3 b k 2 1 + /3/i 2 



(14) 



Note that rcd is our estimate for the scalar field that 
generates the real-space displacement. Next we relate 
^> red to the redshift-spacc galaxy density field: 



V • a rcd 

v y trivial 



rred / a2 

= -IV 2 + ^)^ d (15) 



red 



ired 



It follows from Eq (|T5|) . that q [^viai can ^ e expressed as: 



a rod 

y trivial 
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red 
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From Eq (fl6|) . we can relate q trivial *° 9 trivial as follows 



-.red 

y trivial, z 
-red 

y trivial, x,y 



= (1 



' P)<1 trivial, 



(17) 



- n real 
y trivial, 



where x and y are the directions orthogonal to the line-of- 
sight. Finally, we note that the redshift-space displace- 
ment field is distorted along the z-direction and hence 
we must modulate the field by a factor of 1 + / in the 
z-direction, where / = d\nD/d\na is the logarithmic 
derivative of the growth function. The redshift-space dis- 
placement field is then given by: 



q T z ed : 



l + / 7 red 

^ _|_ n " trivial; z ) 



a rcd = a rcd 

y x,y y trivial; x,y 



(18) 
(19) 



Thus, in redshift space, we displace our data points and 
random points by the displacement field q rcd given by 
Eq (p~9|) . In later sections, we demonstrate that this 
simple reconstruction techniques improves the correla- 
tion between density fields, reduces the measured shift 
in the acoustic scale and reduces the scatter around the 
measured shift. We plan to explore the effects of mov- 
ing the data and random points by different amounts on 
reconstruction in future work. 

3. INVESTIGATING THE CORRELATION OF 
DENSITY FIELDS 

3.1. Introduction to the Propagator 

In linear perturbation theory, the Fourier modes grow 
independently of each other, e.g., there is no mode- 
coupling. In this section, we test how close we are to lin- 
ear theory on a mode-by-mode basis. We look at the cor- 
relation of initial and final (z = 1) density field for each 
mode via the propagator, which compares the amount of 
the initial density field retained in each mode of the final 
density field. We derive the propagator in the presence 
of galaxy bias and com pare it to t he pr opagator for the 
mass case presented bv lSeo et al.l (|2009f ) to estimate the 
additional BAO signal erased by galaxy bias. 



G(k) = 



b x P„ 



a i<(! 



(k) 



'initial 



(*) 



S *final( k ) D ( Z 



oo) 



D{z = \){\+Pn*) 



(20) 



Here 5 initial is the linear theory density field at z = oo, 
5 final is the density field at z = 1, D(z) is the growth 
function, Pi n itiai{k) is the power spectrum of the initial 
density field, and b is the clustering bias for a given biased 
tracer. 

Matter growth is linear on large scales and hence we 
expect that the initial and final density fields should be 
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extremely correlated regardless of the mass tracer used. 
In the absence of non-linear effects and galaxy bias, we 
expect perfect correlation between low redshift and high 
rcdshift density fields which is represented by G(k) = I. 
Figure [2] shows the propagator, G(k) for real and red- 
shift space for HODla. From Figure [21 we note that at 
k w 0.1 h Mpc -1 , the wavenumber around which most of 
the acoustic information is contained, both the real and 
redshift space density fields are correlated extremely well 
with the initial density field: G(k) > 0.9. As expected, 
the non-linear effects and redshift distortions become in- 
creasingly dominant at smaller scales causing the steep 
drop-off in the propagator. We also see that the real 
space density field is correlated with the initial density 
field up to to much smaller scales (larger k) than redshift 
space. The correlation in redshift space is lower on large 
scales than real space due to redshift distortions. 

3.2. Effect of Galaxy Bias on the Propagator 

We next use the HOD models from § 12.21 to evaluate 
impact of galaxy bias on the recovery of the initial density 
field. By using the fact that the initial and final densities 
must be correlated with G(k) = 1 at small fc, we can 
derive the bias for every HOD. In this paper, we use the 
real space G(k) to derive the bias. The bias is reflected 
in the propagator as a multiplicative factor. If we force 
b = 1 in Eq. (f2"U)) , then G(k) will converge to the bias 
at small wavenumber. We fit a constant function to the 
real space propagator, G(k), for k < 0.08 h Mpc -1 to 
compute the bias for each HOD. 

For each HOD, we want to estimate the amount of 
noise increase due to bias relevant to the signal. We de- 
fine an effective number density (n e ff) by Equation (|22[) 
and take any deviation from linear bias into account 
rather than merely using the shot noise. By definition, 
this also parametrizes the degree of the scale dependence 
of the bias. 



=PnOD(k = 0.2 /iMpc" 1 ) (21) 

n eff 

-bias 2 x P Un (k = 0.2 h Mpc" 1 ) 

where Pn n is the linear power spectrum and Phod is the 
measured HOD power spectrum. Table [TJ gives the def- 
initions, HOD parameters, biases, effective number den- 
sity and n eff P(k = 0.2h Mpc -1 ) for each HOD. From 
Table [TJ we see that we span a wide range of effective 
number densities (10~ 3 — 10~ 41 ) and a range of n e jjP 
(2.6-0.6). 

Since we are looking at effects of various HOD models, 
we calibrate our results against the mass case (no HOD), 
which along with three of the twelve HOD models are 
shown in figure [3] for real space (upper panel) and red- 
shift space (lower panel). The three HOD models shown 
in figure [3] represent the spread in the propagators for all 
HOD models and the other HOD models give very simi- 
lar results. In real space, we see the G{k) plots for var- 
ious HOD models are very similar with each other. The 
propagators for different HOD models diverge from each 
other at large k. We attribute this to shot noise which 
depends on the number density of each HOD. Hence, the 
effect of bias on the propagator in real space is very small 



to the acoustic wavenumbers. In redshift space, the sep- 
aration between the HOD models is more pronounced 
because of the Finger of God (FoG) distortions in red- 
shift space. Since different HOD models pick out halos 
at different mass thresholds, the FoG distortions are dif- 
ferent for each HOD model. However, we note that in 
both real and redshift space, the separation between the 
HOD models at the acoustic wavenumbers is small. 

Figure[3]also shows the important result that the prop- 
agators for the HOD models are very similar to the mass 
case. In fact, we see that the biased tracers are slightly 
better than the mas s case. This is con sistent with a sim- 
ilar result found by iNoh et all (|2009f ). We can explain 
this effect as follows. The biased cases replace the dark 
matter halo with a central galaxy and fewer and satel- 
lite galaxies. The central galaxies have the same proper 
velocity as the center of mass of their halo while the satel- 
lite galaxies have the proper velocity of the correspond- 
ing dark matter particles. Therefore, the biased cases 
have experienced less random motion than the mass case. 
Hence, the halos in the mass case are more extended than 
in the biased case. In the mass case, we smear out the 
density field and erase some of the correlation with the 
initial field compared to the biased cases. We see that 
this effect is amplified in redshift space due to redshift 
distortions. While the propagator results are consistent 
with our explaination, w e do not rule out effe cts of halo 
clustering as explored bv lAngulo et al.l ([2008D . 

3.3. Effect of Reconstruction on the Propagator 

We next investigate how reconstruction affects the den- 
sity correlations in both real and redshift space. Figure^] 
shows the comparison of the propagator with and with- 
out reconstruction for 3 different HOD models and the 
mass case. We see that reconstruction restores the corre- 
lation between the initial and final (z = 1) density fields 
to larger wavenumber or smaller scales even for biased 
tracers with a realistic level of shot noise. In fact for 
G{k) = 0.8 the wavenumber increases by a factor of 1.5 - 
2 in real space and redshift space. Reconstruction further 
spreads out the HOD models, which is expected as the 
HOD models have different shot noise levels. As recon- 
struction restores the correlation with the initial density 
fields, we expect that reconstruction will help reduce the 
shifts in the BAO scale and reduce the scatter on mea- 
sured shifts. We discuss these effects of reconstruction in 
the following section. 

4. EFFECT OF GALAXY BIAS ON THE 
ACOUSTIC SCALE MEASUREMENT 

4.1. The Acoustic Scale for Biased Tracers 

In a measurement from a redshift survey, a would be 
the ratio of the acoustic scale in the fiducial cosmology 
to the measured scale, a = 1 means that the measured 
acoustic scale is the same as that predicted by linear 
perturbation theory. In this section, we investigate the 
effects of realistic galaxy bias on the mean a and the 
scatter in a using the method described in § 12.41 The a 
calculated for each HOD and the mass case are given in 
Table [2J Figure [5] plots the scatter in a as a function of 
bias where 6 = 1 case is the mass case. Figure [U plots 
the scatter in a as a function of n e ffP whose values are 
given in Table [TJ 
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For the mass case, iSeo et ail (|2009T > used lower mass 
resolution but higher volume simulations (set G576 in 
that paper) to compute the shift in the acoustic scale 
and found a significant shift at z = 1: a maS s — 1(%) = 
0.13 ±0.04 in real space and 0.16 ±0.06 in redshift space. 
However, the G576 simulations do not have the mass 
resolu tion to sufficiently implement our HODs. ISeo et al.l 
(I200H use another set of simulations with higher mass 
resolution (set G1024). In the G1024 set of simulations, 
they find a smaller, but statistically consistent shift, with 
larger errors: a mass — 1(%) = 0.11 ± 0.16 for real space 
and 0.002 ± 0.23 for redshift space. In this paper, we 
will c ompare ou r resul ts to the G1024 set of simulations 
from ISeo et all (|2009f) to be consistent with the mass 
case. From Tabled in real space, the lower biased HOD 
models (b < 2.6, sets 1,2 and 3) are consistent with the 
linear acoustic scale (a = 1) within the scatter set by the 
sample variance of our simulations. Our highest mass 
HOD models (set HOD4's) shows deviation from a = 1: 
Q:hod4c — 1(%) = 0.57 ± 0.32. We see a similar effect in 
redshift space, where HOD sets 1 and 2 are consistent 
with linear theory (a = 1) and HOD sets 3 and 4 show 
more deviation from linear theory: anoD4c — 1(%) = 
0.79 ± 0.40. We also note that for the biased tracers, 
the scatter in a increases with increasing bias, as we 
are sampling rarer halos and thus shot noise becomes 
more dominant in the highly biased cases. The scatter in 
real space ranges from 0.16% to 0.26% while in redshift 
space, it ranges from 0.22% to 0.43%. Due to redshift 
distortions, we see that the scatter in redshift space is 
greater than in real space. 

However, once we apply reconstruction, we see that all 
the a values are consistent with unity: anoD4c — 1(%) = 
-0.07 ± 0.19 in real space and -0.11 ± 0.20 in redshift 
space. Another important difference between before and 
after reconstruction is the scatter in a. We see that 
in all cases, the scatter in a is reduced by a factor of 
1.5 — 2 with reconstruction than without reconstruction 
in both real and redshift space which is consistent with 
the propagators shown in § 13.31 This confirms that even 
for biased tracers our reconstruction technique is restor- 
ing the information and correlation between initial and 
final density field on the acoustic scale. This is because 
reconstruction corrects for the degradation of the BAO 
signal due to large scale bulk flows and the mass and bi- 
ased tracers bulk flows only differ on small scales. From 
Figure we see that the scatter around a is signifi- 
cantly reduced in real space and in particular redshift 
space once we apply our simple reconstruction scheme. 
We also note that the redshift space with reconstruction 
errors are very close to the real space with reconstruction 
errors. We see this effect for all HOD models and hence 
all biased tracers. This shows us that the reconstruction 
scheme is bias-independent to within the scatter set by 
the sample variance of our simulations. We see a simi- 
lar effect in Fig. |H] where the scatter decreases as n c ff P 
increases and reconstruction reduces the scatter in a. 

4.2. Comparing ajjOD end a mass 

To better understand the effects of galaxy population 
bias on the acoustic scale, we compare the shift in the 
acoustic scale (a—1) to the shift from the mass case. We 
make use of «hod ~ 1(%) vs amass — 1(%) plots. With 
such plots, not only can we analyze how the shifts vary 



between each HOD, but we can also see how they cor- 
relate with the mass case. In these plots, we plot every 
a — 1 value calculated from the resampling method de- 
scribed above. We show one such plot of HOD 2a vs 
mass in Figure [7] From the figure, we see that for both 
real space and redshift space, there is a correlation be- 
tween the shifts measured with the HOD models and 
with the mass case. As expected, we see that redshift 
space has more scatter than real space. However, we see 
that reconstruction not only constrains the distribution 
and thereby reducing the scatter around the mean mea- 
sured shift, but it also increases the correlation between 
the HOD models and mass case. 

4.3. Difference in Acoustic Scale: «hod — o: mass 

In the previous section, we showed that there is a 
strong correlation between the shifts measured from 
the biased tracers and those measured from the mass 
case. This means that we can measure ohod — a mass 
more precisely than chrod alone. This is useful be- 
cause we can measure a mass more pre cisely using sim - 
ulations with poorer mass resolution. ISeo et al.l (|2009l ) 
measured the acoustic scale at z = 1 in redshift space to 
be a mass — 1(%) = 0.158 ± 0.061 before reconstruction 
and a mass — 1(%) = 0.002 ± 0.030 after reconstruction 
using the lower resolution, larger volume simulation set 
(G576). Combining the sets of simulations and precisely 
measuring a mass and anoD — « ma ss allows us to deter- 
mine ckhod- 

The q;hod — awss values are given in Table |3] and 
shown in Figure |SJ From the table and the figure, we 
see some deviation from zero difference for the high 
biased cases which corresponds to low number density 
HOD models for both real and redshift space. The dif- 
ference in redshift space is slightly larger than in real 
space: «hod4c — awss(%) = 0.69 ± 0.24 in real space 
and 0.79 ± 0.31 in redshift space. However, the differ- 
ence between the shifts from biased tracers and the mass 
case, even for the highest biased cases, is about 2.7c. We 
also note that the scatter around the difference increases 
with increasing bias. 

When we apply our reconstruction scheme, we see that 
there is no difference between the shifts measured by bi- 
ased tracers and the mass case for both real and redshift 
space: «hod4c ^ awss(%) = —0.03 ± 0.16 in real space 
and -0.05 ± 0.16 in redshift space for HOD4c. The dif- 
ference between the shifts is consistent with zero with 
errors < 0.08% for HOD models 1 and 2. Equally impor- 
tant, we also note that reconstruction reduces the scatter 
around the difference in the shifts by a factor of 1.5 — 2 
depending on the HOD model. Also, the scatter around 
the difference in the shifts after reconstruction are the 
same for both real and redshift space. This is consistent 
with our previous result that reconstruction accounts for 
non linear structure formation and these flows for the 
biased cases and mass tracers differ only at small scales. 

4.4. Slopes of (Jhod — 1 vs a mass — 1 Distribution 

From Figure [71 we see that the qhod and a maS s values 
are correlated. In this section, we investigate the slope 
of the «hod — 1 vs a mass — 1 distribution. If the HOD 
models and mass cases are perfectly correlated, then we 
expect the underlying slope to be unity. However, there 
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are two effects than can lead to non-unity slopes. First, 
if the shifts are slightly uncorrelated between the mass 
and HOD models, we can see a non-unity slope. Sec- 
ond, even if the shifts are perfectly correlated (underly- 
ing unity slope), a difference in the scatter around the 
mean between the mass case and the HOD models can 
lead to a non-unity slope. 

In Figure [3 we compute the slope by using a linear 
least squares fits. However, we know that this slope and 
the error for the slope are not correct since the linear least 
squares method assumes that all the scatter is around the 
y-axis variable, «hod — 1- We test the statistical signif- 
icance of a non-unity slope by running multiple Monte- 
Carlo simulations with an underlying distribution with 
unity slope. We then compare the results of the Monte- 
Carlo simulations to our measurements. In detail, we use 
Monte Carlo simulations given by Eq. (f2"2"j) , which has an 
underlying slope of unity to calculate the measured slope. 
Our model is: 

a maS s,MC - 1 = Gaussian(/x = 0, a = cr mass ), 

«HOD,MC - 1 = "mass.MC + GauSSian(^i = 0, a = CTHOD-n 

(22) 

where <r mass is the scatter around a mass from Table [5] 
and <7HOD-ma SS is the scatter around a H OD-mass horn 
Table [3J We use the same resampling method as used 
for our simulations to generate 1000 subsamples by av- 
eraging over 22 randomly selected values from 44 values 
for each subsample given by Eq. (f22|). In order to com- 
pute the mean Monte Carlo slope and the scatter around 
this slope, we create 500 realizations of 1000 subsamples 
and use linear least squares fit. In Table 21 we give the 
slopes for «hod vs a mass distribution for a mass at z = 1 
and z = oo for four different HOD models with the other 
HOD models giving very similar results. We see that for 
both redshift space and real space, the low bias and high 
number density HOD models (HODs 1 & 2) are all con- 
sistent with an underlying slope of unity. At the high 
bias and low number density HOD models (HODs 3 & 
4) , we start to see some deviation from unity slope at the 
1.5a — 2.0cr level for z — 1 mass values. However, with 
reconstruction, we see that all HOD models are consis- 
tent with an underlying unity slope for both z = 1 and 
z = oo mass values within ler. 

4.5. Fisher Matrix Scatter Estimates 



ative to the Fisher matrix code from ISeo fc Eisenstein 
' 20071) as shown by iTakahashi et al.l (|2009t) fc ISeo et al. 
20Q9I). This is because in our A^-body analysis we fit 



ISeo & Eisenstein! (|2003ll2007l ) used Fisher matrix anal- 
yses to predict the scatter in the acoustic scale available 
in surveys of a given number density and bias. These 
models depend on certain assumptions, such as Gaus- 
sianity of the density field up to a cut-off wavenum- 
ber, that we can check with A^-body simulations. Here 
we compare the scatter in the acoustic scale for our 
HOD models to those predicted by the Fisher matrix 
model. The Fisher matrix code takes in the number 
density (Eq (|22jl). a 8 (Table [I]), Ej_, Ey and /3. For 
our cosmology, at z = 1, Ej_ = 5.26 is the rms ra- 
dial displacement across the line of sight and Eii = 5.26 
(real space) and 9.55 (redshif t space) is the rm s displace- 
ment along the line of sight (|Seo et al.ll2009D . When we 
apply reconstruction, we choose T,± t p> econ = Ej_/2 and 
^\\,Recon = We know that the iV-body analysis 

tends to overestimate the scatter in redshift spaec rel- 



a spherically averaged power spectrum, we do not opti- 
mally extract the two-dimensional information. Hence, 
we expect the scatter around the acoustic scale from the 
A^-body simulations to be a factor of vL16 = 1.08 greater 
than the Fisher matrix estimates at z = 1. 

Table [2] provides the a — 1 results from our simulations 
and the Fisher matrix estimates for the scatter are given 
in brackets. We see that for the low biased HOD models 
(HODs 1 and 2), the scatter around the mean measured a 
is very close to the Fisher matrix estimates for real space 
and redshift space. For the high biased HOD models 
(HODs 3 and 4), we see that the measured scatter is 
larger than the the Fisher matrix estimates by 20% — 
30%, which is more than expected. It is unclear if the 
HOD sets 3 and 4 show larger scatter due to high bias or 
low n e ffP. However, once we apply reconstruction, we 
have very good agreement between the measured scatter 
and the Fisher matrix estimates for all HOD models. 

) 

4.6. Comparing a rca i and a re dshift 

Since we measure galaxies in redshift space, we do not 
know the real-space positions or velocities of galaxies. In 
this section we see how well the shifts measured from 
redshift space are correlated with real space values. We 
use the linear least squares fit to c omp ute the slopes 
of the distribution as done in § 14.21 and 14.41 However, 
we know that not only are the shifts measured in real 
and redshift space correlated, but the scatter around the 
shifts are also correlated. In Table [5l we quote the lin- 
ear least squares slopes for a re dshift vs a roa i distribution 
with and without reconstruction. We expect a non-unity 
slope from for the distributions since we know the scat- 
ter around the rcdshift-spacc values arc larger than the 
real space values. We see that this is true from the val- 
ues in Table [3 We also see that reconstruction decreases 
the slope towards unity. We expect this behavior since 
reconstruction reduces the noise in redshift space that 
is uncorrelated with real space and vice- versa. Figured] 
shows the a re ai vs o; re dshift plot for HOD2a and the mass 
case. From such plots, we can compare the shifts derived 
from real and redshift space. We see that reconstruction 
increases the correlation between real and redshift space 
and reduces the scatter for both sets of values. Wc also 
see that the HOD and mass cases are very similar in their 
real and redshift space correlation after reconstruction. 
The other HOD models give very similar results as shown 
in Table 

4.7. Effects of Variations in Bias and Reconstruction 
Smoothing Scale 

When wc implement our reconstruction scheme out- 
lined in § 12.51 we use a smoothing scale of R = 
XAhT 1 Mpc. We change the smoothing scale to R = 
20/) _1 Mpc and rerun our reconstruction scheme. We see 
no change in either a or a a when we change the recon- 
struction smoothing scale. Thus, we see that reconstruc- 
tion gives consistent answers over a range of smoothing 
scales. We also do not see any differences in «hod — admass 
values when using the larger smoothing scale for «hod ■ 

To implement reconstruction, we need to estimate the 
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value of the galaxy bias for a survey. We explore the ef- 
fects of inputing a wrong bias for calculating the acoustic 
scale and reconstruction techniques. We use biases that 
are incorrect by ±10%, ±20% and ±40%. The results are 
plotted in Figure ITOl From the figure, we sec that chang- 
ing the bias obviously does not affect the propagator be- 
fore reconstruction. After reconstruction, the propagator 
is very close to the actual bias case at the acoustic scale 
for ±10%. We see that the ±40% case shows moderate 
deviation from the correct bias case. Thus, we conclude 
that our reconstruction scheme is insensitive to a large 
range (up to ±30%) of incorrect input bias. From future 
weak lensing and cluster surveys, we expect to know the 
amplitude of matter clustering well enough for us to know 
the bias to less than 10%. Hence, we expect this effect 
to be negligible. 

5. CONCLUSION 

We use high force resolution iV-body simulations with 
a total volume of 44/i -3 Gpc 3 to analyze the effects of 
galaxy bias on the Baryon Acoustic Oscillation (BAO) 
scale. In particular, we apply a variety of Halo Occu- 
pation Distributions (HODs) to simulate galaxy popula- 
tions and galaxy bias for upcoming dark energy surveys. 
We use a variety of techniques to analyze the shift in 
the acoustic scale for different galaxy populations and 
co mpare the derive d shifts with the mass case presented 
in iSeo et all (|2009f ) . We ext end the simple reconstr uc- 
tion technique introduced by lEisenstein et all (|2007t ) to 
biased tracers and use it to remove effects of large scale 
gravitational flows and thereby preserving the linear the- 
ory density field. 

We use the propagator to look at the damping of the 
BAO signal due to non-linearities for different HOD mod- 
els. We analyze our results in both real space and redshift 
space and apply reconstruction to both. Wc find that the 
damping is more in redshift space compared to real space. 
We find that galaxy bias has little effect on the propa- 
gator at BAO wavenumbcrs. In fact, the propagator for 
even the highest biased cases are close to the mass case. 
Reconstruction restores correlation between the initial 
and final density fields for both real space and redshift 
space. Hence, reconstruction removes most of the degra- 
dation due to non linear structure formation. We find 
that with reconstruction G(k) w 0.9 for k < 0.3h Mpc -1 
in real space and k < 0.2h Mpc^ 1 in redshift space at 
z = 1. The spread in HOD models after reconstruction 
can be accounted by the different levels of shot noise and 
redshift distortions for each HOD. 

We use our power spectrum fitting to measure the 
shift in the acoustic scale. We detect a mild shift 
in the acoustic scale in real space for our high bias 
(b > 3) HOD models. For the most biased HOD model, 
o;hod4c — 1(%) = 0.57 ± 0.32. In redshift space, the 
low biased cases are consistent with no shift at la while 
the high biased cases show a shift less than the 2a level: 
«HOD4c — 1(%) = 0.79 ± 0.40. However, once we use 
reconstruction, we do not detect any shift for any HOD 
in either real space or redshift space. Also, we see that 
reconstruction reduces the scatter around the measured 
shift by a factor of 1.5 — 2. This confirms our results from 
the propagator that reconstruction undoes the degrada- 
tion to the BAO signal caused by large scale bulk flows. 

In addition to the shift in the acoustic scale, we also 



compare the shift to the mass case. We look at the 
difference in the shift computed with the HOD mod- 
els and the mass case and look at the scatter around 
this difference. We see that for the low biased HOD 
models, wc sec no difference between «hod and a mass 
within la: aHOD2b — a mass (%) = —0.05 ± 0.09 in real 
space and 0.05 ± 0.11 in redshift space. The high bi- 
ased HOD models show a difference at the 1.5 — 2.7a 
level: «hod4c — Q? mass (%) = 0.69 ± 0.25 in real space and 
0.79 ± 0.31 in redshift space. Once we apply reconstruc- 
tion, the difference between the shifts are consistent with 
within la. Reconstruction also reduces the errors on 
the difference of the shifts by a factor of 1.5 — 2 depending 
on the HOD model in real and redshift space. In sum- 
mary, the acoustic scale in redshift space matches that 
of linear theory to within 0.08% for HOD models 1 & 2, 
0.12% for HOD models 3 and 0.16% for HOD models 4. 

We also test the correlation between the shifts mea- 
sured with HODs and without HODs using ohod vs 
admass plots and the slopes of these distributions. From 
both these analyses we conclude that there is strong cor- 
relation between the shifts. We also show that recon- 
struction not only reduces the scatter, but also improves 
the correlation. With reconstruction, our results between 
shifts from HOD models and the mass case consistent 
with an underlying slope of unity within error bars. Wc 
use a similar method to analyze the correlation between 
shifts from real space and redshift space. We find that, 
as expected, reconstruction strengthens the correlation 
between real space and redshift space by accounting for 
redshift distortions. 

We use the Fisher matrix code by ISeo fc Eisensteinl 
(|2007| ) to compare the measured scatter around the 
acoustic shift to theoretical predictions. We find that for 
the low biased HOD models (HODs 1 and 2), we have 
very good agreement between the measured scatter and 
the predicted scatter. For the high biased HOD mod- 
els (HODs 3 and 4) however, we find that the measured 
scatter is larger than the predicted scatter by 20% — 30%. 
However, after reconstruction, we have good agreement 
between the measured and predicted scatter for all HOD 
models. 

Finally, we look at effects of varying the reconstruction 
smoothing scale and input galaxy bias. In our analysis, 
we have used a smoothing scale of 14/i _1 Mpc. We in- 
crease that to 20h~ 1 Mpc and we do not find any signifi- 
cant differences in the acoustic shift or the scatter around 
the shift. We considered the impact of misestimation of 
the bias parameter in the reconstruction method. Vary- 
ing by 10%, 20%, and 40%, we found mild degradations 
in the propagator only for the largest case. We expect 
future surveys to know the bias to better than 10%, more 
than adequate for reconstruction. 

In conclusion, we investigate the effect of galaxy bias 
on the acoustic scale. We find that the effect of galaxy 
bias in redshift space relative to the mass case without re- 
construction is about 0.1% for the low bias cases (b < 3), 
growing to 0.3% in the most extreme case. With our sim- 
ple reconstruction scheme, this effect is consistent with 
no shift with errors less than 0.08% for the low biased 
cases and less than 0.16% for the high biased cases. Cur- 
rent surveys such as the Sloan Digital Sky Survey (SDSS) 
III baryon oscillations spectroscopic survey (BOSS), the 
WigglcZ dark energy survey, and the Hobby-Ebcrly tele- 
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scope dark energy experiment (HETDEX) will measure 
the acoustic peak to about 1%, 2.5%, and 2% precision 
respectively. While our simulations probe the effects 
of galaxy bias at higher redshift than current surveys, 
our results suggest that this effect will not be noticeable 
in these surveys. However, future surveys such as Big- 
BOSS, Joint Dark Energy Mission (JDEM), and Euclid 
will measure the acoustic scale to a precision close to 
the cosmic variance lim it of about 0.1% out to z = 2 
(|Seo fc Eis eustein 2007). Our current error bars are ap- 



proaching these levels, but we plan to running more high 
resolution simulations to lower redshifts with larger vol- 
umes in order to reduce our error bars. Future work will 
also look at effects of cut-sky and non-periodic boundary 
conditions on reconstruction. 

We thank Martin White for helpful discussions. KM, 
DJE, JE, and XX are supported by NSF AST-0707725 
and by NASA BEFS NNX07AH11G. HJS is supported 
by the DOE at Fermilab. 
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Fig. 1. — The average power spectrum from 44 simulations (HOD la, real space) divide d by the no wiggle power spectrum, P(k)/ Pnwl(k). 
The no-wiggle power spectrum represents the power spectrum without the BAO peaks Eiscnstcin & Hu (1998). Thus, we clearly see the 
BAO peaks appear in the power spectrum. We also start to see the increasing power on small scales (large k) as non-linear growth starts 
to dominate. Plotted in the dashed line is the average power spectrum after reconstruction. We see that our reconstruction scheme reduces 
the effect of non-linear structure formation on small scales and restores information into the acoustic peaks. The dot-dashed line rcprcsentcs 
the linear power spectrum divided by the no-wiggle power spectrum. 




Fig. 2. — The correlation between the initial and final density field as described by Eq. 1201 for HOD la. The solid line represents real 
space while the dashed line represents redshift space. We see that redshift distortions reduce the correlation between the initial and final 
(z = 1) density fields. We are most interested in the correlation between the density fields at BAO wavenumbers (k ks 0.09 — 0.2h Mpc -1 ). 
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Fig. 3. — The propagator, G(k) for three HOD models and the mass case. The top panel shows real space and the bottom panel shows 
rcdshift space. We see that while redshift distortions reduce the correlation between the density fields, they do not spread the individual 
HOD models from the mass case. 
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k [h/Mpc] 




k [h/Mpc] 

Fig. 4. — The propagator G(k) for three HOD models and includes the effects of reconstruction. The top panel shows real space and 
the bottom panel shows redshift space. In both panels, the lines maintaining higher correlation at smaller scales are for the cases after 
reconstruction. We see that in both real and redshift space, reconstruction restores information on smaller scales. 



14 



8 

a 

•FN 

u 

QJ 
OS 

in 



0.5 
0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 



L 
0.5 



A real space 
□ redshift space 
A real + recon 
■ redshift + recon 

mass 



galaxies 



A 
1 



Ob 

^A 



2x10 3 1x10 



-3 



1.5 



2 

bias 



□ 
□ □ 



2>^ 







-4 



1x10 



-4 



2.5 



3.5 



Fig. 5. — Scatter in the a values as a function of bias, bias = 1 refers to the mass case. The plot also shows the effective number density 
( n eff) m units of h 3 Mpc - 3 for each HOD set. From the plot, we see that reconstruction dramatically reduces the error in a for all HOD 
models in both real and redshift space. We note that, the redshift space with reconstruction errors are very close to the real space with 
reconstruction errors even for the highest biased case. This shows that reconstruction works well for all biased tracers. 
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Fig. 6. — Scatter in the a values as a function of n e jjP. From the plot, we see that the scatter in a reduces as we go to higher n e jP. 
Also, we see that reconstruction reduces the scatter in alpha consistent with the results we found and with Fig. \5\ 
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Fig. 7. — ctHOD vs amass plot for HOD 2a for real and redshift space, with and without reconstruction. From such a — a plots, we can 
compare the a values from the HOD models to the mass case. The crosses represent the mean a — I and the scatter around the mean. 
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Fig. 8. — ohod — ctmass values for various HOD models corresponding to different biases and different number densities. We see some 
evidence for a shift between the most biased HODs and the mass case. However, we do not see this effect after we apply reconstruction. 
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Fig. 9. — «rcal vs "^redshift plot for HOD2a and the mass case, with and without reconstruction. From such a — a plots, we can study 
the correlation between the shifts derived from real and redshift space. The crosses represent the mean a — 1 and the scatter around the 
mean a. 



19 



10% 


less bias 


10% 


more bias 


20% 


less bias 


-20% 


more bias 


40% 


less bias 


40% 


more bias 


Correct bias 




10 



-1 



k [h/Mpc] 



10* 



Fig. 10. — The effect of inputing a wrong bias in the reconstruction scheme. The acoustic scale is not affected upto a 30% change in bias. 
We see moderate deviations from the correct bias case for the 40% change in bias case. 



20 



TABLE 1 

Halo Occupation Distributions (HODs) and Their Properties. 



Model Total Number Satellite M cen M aat bias n cff (h 3 Mpc -3 ) n cff P(k = 0.2 h Mpc -1 ) 
of Galaxies Fraction (%) 



HODla 


2 


X 


10 6 





1.37 




1.75 


1.14 


X 


10" 3 


2.59 


HODlb 


2 


X 


10 6 


5 


1.37 


0.92 


1.79 


1.10 


X 


10~ 3 


2.59 


HODlc 


2 


X 


10 6 


10 


1.37 


0.44 


1.82 


1.05 


X 


10" 3 


2.58 


HOD2a 


1 


X 


10 6 





2.52 




1.98 


6.79 


X 


10" 4 


1.97 


HOD2b 


1 


X 


10 6 


5 


2.59 


1.52 


2.03 


6.37 


X 


10~ 4 


1.93 


HOD2c 


1 


X 


10 6 


10 


2.74 


0.74 


2.08 


5.80 


X 


10" 4 


1.86 


HOD3a 


3 


X 


10 s 





6.19 




2.46 


2.67 


X 


10~ 4 


1.20 


HOD3b 


3 


X 


10 s 


5 


6.41 


3.06 


2.52 


2.45 


X 


10~ 4 


1.15 


HOD3c 


3 


X 


10 s 


10 


6.70 


1.48 


2.58 


2.23 


X 


10~ 4 


1.09 


HOD4a 


1 


X 


10 s 





14.40 




3.17 


9.10 


X 


10~ 5 


0.68 


HOD4b 


1 


X 


10 s 


5 


14.80 


5.98 


3.23 


8.22 


X 


10" 5 


0.63 


HOD4c 


1 


X 


10 5 


10 


15.48 


2.88 


3.29 


7.51 


X 


10" 5 


0.60 


a M ce „ is 


given 


in 


units of 10 12 


h- 1 


M© and M aa t 


in units of 10 14 


H^Mq 







b O-bias = "mass * D ( Z = l )/ D ( z = °) with ""mass = 0.817. 

c n e fF is given by Eq. (122t . 



TABLE 2 

a — 1 (%) values. 



Model 


Real Space 


Real Space 


Redshift Space 


Redshift Space 




(Theoretical) 


with Reconstruction 


(Theoretical) 


with Reconstruction 



mass 


-0 


11 


± 





16 


(0 


15) 


-0 


00 


±0 


10 


(0 


09) 


-0 


01 ±0 


23 


(0 


21) 


-0 


06 ±0 


11 


(0 


10) 


HODla 


-0 


10 


± 





16 


(0 


16) 


-0 


01 


±0 


12 


(0 


12) 


-0 


05 ±0 


22 


(0 


22) 


-0 


03 ±0 


13 


(0 


11) 


HODlb 


-0 


08 


± 





17 


(0 


16) 


+0 


00 


±0 


12 


(0 


12) 


+0 


10 ±0 


22 


(0 


22) 


+0 


00 ±0 


13 


(0 


14) 


HODlc 


-0 


11 


± 





17 


(0 


17) 


+0 


05 


±0 


12 


(0 


12) 


+0 


04 ±0 


22 


(0 


22) 


-0 


05 ±0 


13 


(0 


14) 


HOD2a 


-0 


09 


± 





10 


(0 


17) 


-0 


03 


±0 


12 


(0 


13) 


+0 


10 ±0 


24 


(0 


23) 


-0 


05 ±0 


13 


(0 


15) 


HOD2b 


-0 


07 


± 





10 


(0 


18) 


-0 


05 


±0 


12 


(0 


13) 


+0 


14 ±0 


24 


(0 


24) 


-0 


06 ±0 


13 


(0 


15) 


HOD2c 


-0 


16 


± 





10 


(0 


18) 


-0 


11 


±0 


12 


(0 


13) 


+0 


04 ±0 


23 


(0 


24) 


-0 


07 ±0 


13 


(0 


15) 


HOD3a 


-0 


06 


± 





27 


(0 


20) 


-0 


11 


±0 


17 


(0 


15) 


+0 


26 ±0 


33 


(0 


27) 


-0 


14 ±0 


18 


(0 


17) 


HOD3b 


+0 


07 


± 





26 


(° 


20) 


-0 


07 


±0 


17 


(0 


15) 


+0 


39 ±0 


32 


(0 


27) 


-0 


07 ±0 


18 


(° 


16) 


HOD3c 


+0 


06 


± 





27 


(0 


20) 


-0 


11 


±0 


17 


(0 


15) 


+0 


42 ±0 


33 


(0 


27) 


-0 


03 ±0 


18 


(0 


16) 


HOD4a 


+0 


12 


± 





33 


(o 


25) 


-0 


15 


±0 


20 


(0 


19) 


+0 


62 ±0 


40 


(0 


33) 


-0 


15 ±0 


21 


(o 


22) 


HOD4b 


+0 


11 


± 





33 


(0 


25) 


-0 


27 


±0 


20 


(0 


20) 


+0 


69 ±0 


42 


(0 


33) 


-0 


23 ±0 


21 


(0 


23) 


HOD4c 


+0 


57 


± 





32 


(0 


26) 


-0 


07 


±0 


19 


(0 


20) 


+0 


79 ±0 


40 


(0 


31) 


-0 


11±0 


20 


(0 


23) 



a Mass values taken from lSeo et al.l |2009). 

b The real and redsh ift space values given in the bracket are the theoretical Fisher matrix predictions 
for the scatter in a l lSeo & Eisenstcin 20071 ) 



TABLE 3 

"HOD - Ctitass (%) VALUES. 



Model 


Real Space 


Real Space 


Redshift Space 


Redshift Space 






with Reconstruction 




with Reconstruction 


HODla 


+0.02 ±0.08 


+0.08 + 0.05 


+0.06 ± 0.09 


+0.03 + 0.06 


HODlb 


+0.00 ±0.07 


+0.05 ± 0.05 


+0.05 ± 0.08 


+0.01 + 0.07 


HODlc 


+0.03 + 0.07 


+0.09 ± 0.06 


+0.10 + 0.08 


+0.06 + 0.07 


HOD2a 


+0.02 ± 0.09 


+0.06 + 0.07 


+0.11 + 0.10 


+0.01 + 0.07 


HOD2b 


-0.05 ± 0.09 


-0.01 + 0.07 


+0.05 + 0.11 


-0.02 + 0.08 


HOD2c 


+0.04 ± 0.09 


+0.04 ± 0.07 


+0.15 + 0.11 


+0.00 + 0.08 


HOD3a 


+0.05 + 0.14 


-0.05 + 0.12 


+0.26 + 0.18 


-0.08 + 0.12 


HOD3b 


+0.18 + 0.16 


-0.01 + 0.11 


+0.43 ± 0.20 


+0.03 + 0.12 


HOD3c 


+0.18 + 0.14 


-0.02 + 0.12 


+0.39 + 0.18 


-0.06 + 0.12 


HOD4a 


+0.53 + 0.24 


-0.06 + 0.16 


+0.63 + 0.29 


-0.09 + 0.16 


HOD4b 


+0.52 + 0.25 


-0.17 + 0.16 


+0.70 + 0.31 


-0.16 + 0.17 


HOD4c 


+0.69 + 0.25 


-0.03 + 0.16 


+0.79 + 0.31 


-0.05 + 0.16 



+ Mass values taken from lSeo et al.l (2009). 



TABLE 4 

Slopes of the qhod vs a m ass distribution. 





Mass values at z = 1 


Mass values at z = oo 


HOD 


Measured Slope 


Monte Carlo Slope Measured Slope 


Monte Carlo Slope 


Heal Space 


HODlb 


0.96 


1.00 ±0.06 


1.07 


1.00 ±0.26 


HOD2b 


1.04 


0.99 ± 0.09 


1.02 


1.00 ±0.31 


HOD3b 


1.42 


1.00 ±0.15 


0.89 


1.02 ± 0.46 


HOD4b 


1.42 


1.00 ±0.24 


0.71 


0.98 ±0.57 






Real Space With Reconstruct: 


ion 




HODlb 


1.11 


0.99 ±0.09 


0.99 


1.00 ±0.15 


HOD2b 


1.07 


1.00 ±0.12 


0.98 


1.00 ±0.17 


HOD3b 


1.31 


1.01 ±0.21 


0.91 


1.02 ±0.27 


HOD4b 


1.28 


1.01 ±0.27 


0.98 


0.99 ± 0.34 


Rcdshift Space 


HODlb 


0.88 


1.00 ±0.06 


1.14 


1.03 ±0.35 


HOD2b 


0.91 


1.00 ±0.07 


1.12 


0.98 ±0.39 


HOD3b 


1.15 


1.00 ±0.12 


1.02 


0.94 ±0.57 


HOD4b 


1.22 


1.00 ±0.21 


0.26 


1.05 ± 0.79 


Redshift Space With Reconstruction 


HODlb 


1.14 


1.00 ±0.10 


0.96 


1.02 ±0.19 


HOD2b 


1.08 


1.00 ±0.12 


0.96 


0.99 ±0.19 


HOD3b 


1.31 


1.01 ±0.20 


0.93 


1.00 ±0.28 


HOD4b 


1.28 


1.01 ±0.28 


0.79 


0.98 ± 0.38 



TABLE 5 

Slopes for the a re d B hift vs « re ai distribution. 



Model 


No Reconstruction 


With Reconstruction 


mass 


1.30 


1.00 


HODla 


1.19 


1.08 


HODlb 


1.18 


1.10 


HODlc 


1.21 


1.08 


HOD2a 


1.17 


1.01 


HOD2b 


1.15 


1.04 


HOD2c 


1.13 


1.00 


HOD3a 


1.17 


1.04 


HOD3b 


1.14 


1.05 


HOD3c 


1.16 


1.03 


HOD4a 


1.18 


1.05 


HOD4b 


1.19 


1.01 


HOD4c 


1.20 


1.02 



